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Abstract 

We evaluate perturbative next-to-leading order corrections to the hard exclusive leptoproduction 



of 7r + mesons on a transversely polarized proton target. A model dependent study shows that 
j> \ these corrections can be large. We analyze the scale dependence and explore the Brodsky-Lepage- 

•rH , 

^ . Mackenzie scale setting procedure. Although the predictions for the cross section suffer from 
5^ ■ 

theoretical uncertainties, the transverse nucleon single spin asymmetry turns out to be a rather 

_ _ 1 

stable observable since higher order effects approximately cancel there. 
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1 Exclusive meson production and QCD factorization. 



The hard exclusive leptoproduction of a meson M from a nucleon target N, 

£(k)N(P 1 ) - £'(k')N'(P 2 )M(q 2 ) (1) 

is a promising process to test our understanding of QCD in exclusive reactions, as well as a means 
for studies of the properties of nucleon to hadron transitions, N — > N', with N' being a baryon 
from an SUf(3) multiplet. If the intermediate photon is longitudinally polarized and has a large 
virtuality Q 2 = —qf, the photoproduction amplitude ^ L N — ► N'M factorizes into a convolution 
of three parts [I], see Fig. [I], 

A= [ du f dx 4> ff (u\n)T ff (u,x, £| Q,n)Aff>{x,£, A 2 + . . . , (2) 

f,f'=u,d,s 1 

where the ellipsis stand for power suppressed contributions in 1/Q. The hard subprocess Tfp 

encodes the short distance dynamics of the parton scattering and can be consistently calculated 

in QCD perturbation theory as a series in the strong coupling a s . The other two blocks, <fi 

and A, are universal, i.e. process independent, and embody the long-distance physics. is a 

conventional leading twist meson distribution amplitude. It describes the minimal Fock component 

of the meson wave function with two quarks having momentum fractions u and 1 — u w.r.t. the 

meson momentum, respectively. A is a flavour nondiagonal generalized parton distribution (GPD). 

Its definition, taking apart the flavour transition, differs from the conventional Feynman parton 

densities by the presence of a non-zero momentum flow, A = P 2 — P± = qi — q 2 , in the t-channel. As 

a result, the GPD is a complicated function of three variables, the s-channel momentum fraction 

x, its t-channel counterpart^ rj ~ A + , and A 2 . The appearing phase space picture is rich and 

results into a trinity interpretation of GPDs: they share common properties with forward parton 

densities, distribution amplitudes, and form factors. 

In the present study we will concentrate on the next-to-leading order (NLO) analysis of the 

7r + -pro duct ion from the proton. Thus we set N = p, N' = n and M = ir + in Fig. [l|. This reaction 

is an issue of intensive experimental studies by HERMES and Jefferson Lab in relation to 

the measurement of the helicity-flip GPDs accessible in this process. Our consequent presentation 

is organized as follows. In the next section we introduce the definitions of all basic ingredients 

entering the amplitude (|2]) and calculate the electroproduction cross section for the transversely 

polarized proton. In section |] we present one-loop corrections for the hard-scattering amplitude. 

In section f| we introduce simple models for the GPDs and then we discuss the issue of scale 

1 Note that in the kinematics, we are considering, the skewedness parameter 77 is approximately equal to a 
(negative) generalized Bjorken variable £, to be specified below. 
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Figure 1: Factorization of the meson leptoproduction at — q\ — > oo into hard-scattering amplitude 
T and non-perturbative functions A, the generalized flavour changing parton distribution, and </>, 
the distribution amplitude of the outgoing meson. 

setting. After this we give numerical estimates of the differential cross section and transverse 
single spin asymmetry. Finally, we summarize. 

2 Amplitude and cross section for j L p — > n + n. 

The hadronic part of the 7p — > n + n process is given by the Fourier transform of the matrix element 
of the electromagnetic current J M = ^A7ia em Y,i Qi^il^i 

J d A xe- i ^ x (n + (q 2 )n(P 2 )\J,(x)\p(P 1 )) = t(2n) 4 5^ ( Ql + P 1 - q 2 - P 2 ) A? , (3) 

with qi = k — k' being the difference of incoming and outgoing lepton momenta. Its leading 
term in 1/Q 2 is picked up by the contraction with the longitudinal polarization vector e L . A 
straightforward leading twist calculation gives for the amplitude 

_ i .J -l „ Aud 

A; + = VA^ — ^jJ du dxMu)T ud (u,x,0 1 —t r (x,Z,A 2 ) + O(l/q 3 ), (4) 
p N c q z Jo J-i q ■ P v 7 

where q = | (qx + q 2 ), P = P\ + P 2 and A = q 1 — q 2 = P 2 — P x . The scaling variable is a 

generalized Bjorken variable £ = — ~ — Here in the last equality we neglected the pion 

mass, = 0, and power suppressed corrections in the nucleon mass M 2 /(-q 2 ). The Lorentz 

structure j M = 2q fl + 3£P M , appearing in the twist-two part of A p + , is gauge invariant to the 

twist-four accuracy, i.e. qi^A 1 ^ = \Q+ ^) ^ + ~ 0- 

We introduced in Eq. @ two non-perturbative objects, the pion distribution amplitude 

(n + (q 2 )\u(y) lpl5 d(z)\0) = -iq 2p f % £ du e iu ^ +i ^ u ^ ^(u) , (5) 

with the decay constant f n ~ 132 MeV, as well as A u p d expressed in terms of the off-forward matrix 
element of the non-local flavour-changing operator 

(n(P 2 )|%) 7p 75^)|p(Pi)) = j\xe^ A+xP >y + ^ A - xP > z A u p d (x,^,A 2 ). (6) 
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Figure 2: Hard-scattering coefficient function at leading order. 

It gives rise to flavour non-diagonal GPDs, which enter as form factors 

Af{x,i, A 2 ) = h p H ud (x, £, A 2 ) + e p E ud (x, £, A 2 ), (7) 

in front of the Dirac structures 

/V = 0n(P 2 )7p75^p(Pi), e p = ^C/ ri (P 2 ) 75 C/ p (P 1 ). (8) 

We can safely neglect the mass difference of proton and neutron and set in the following M = 
M p = M n . Using the isospin symmetry of strong interactions, one can reduce the above GPDs to 
the conventional flavour diagonal proton matrix elements |3J 

A ud = A uu _ A dd ^ 

Evaluating the tree diagrams presented in Fig. ^ we get the known result for the hard amplitude 
T ud to leading order (LO) accuracy in the QCD coupling [|, |. [/J, 

T ud {u, x, = C F a s \ | ^ r- - , Qd v 1 + 0(a 2 s ). (10) 

Z \(l-w)(l-f-i0) u(l + f-i0)J 

Here the quark charges are Q u = 2/3 and Qd = —1/3. 

Let us remark on similarities of the present consideration to the one of the pion form factor since 
we will use them below. Contributing diagrams can be decomposed into two sets, where in the first 
(second) one the photon is attached to the (anti-) quark line. The momentum of the initial and final 
(anti-) quark is given in the collinear approximation by pi = — ^^A (p[ = —^-A^j and P2 = uq2 
(p' 2 = (1 — u)q 2 ), respectively. In leading twist approximation both of these sets separately respect 
current conservation. Obviously, if we formally replace — > v (then — > 1—v) with < v < 1, 
we reduce the result ( PD to the kinematics of the pion form factor. Of course, we can also recover 
the amplitude (0) from the pion form factor result by the substitution v — > In addition 

we have to keep track of the emerging imaginary part, which is absent in the pion form factor. 
It develops now in the region \x\ > £ and can be easily restored. Using this analogy, the NLO 
corrections of T ud will be evaluated in section [| 
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Figure 3: The reference frame for the cross section is defined as a target rest-frame and the z-axis 
directed along the three-vector of the virtual photon. The outgoing hadrons, pion and neutron, 
are in the same plane (due to the momentum conservation), which form an angle (p with the 
transverse polarization vector S of the target. 

Now we turn to the calculation of the cross section for electroproduction of pions from the 
proton target £(k)p(Pi) — > £(k')n(P 2 )iT + (q 2 ) . It is given by 



da 7 



1 



L A 77 

A 4 
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dLIPS 3 , (11] 



where the hadronic amplitude (|j) is contracted with the leptonic current (which includes the 
photon propagator) 



i 

Q 

and the three-particle phase space volume reads 



Lfj, = —V^a^uik'^f.uik), (12) 



.LIPS3 = (2^>(* + Pi-k' -P 2 - *) J { i )3 J 2 ^%r d3) 

In the rest-frame of the target with the z-axis chosen along the momentum of virtual photon, as 
shown in Fig. |3], we get the following four-fold cross section 



da* + 1 x B y 2 ( i A M 2 xlY 1/2 



dQ 2 dx B d\A 2 \dip 2(4vr) 4 Q A \ Q 



2 



L A 77 



2 



(14) 



where we use the variables Q 2 = —qf, xb = — <Zi/(2gi ■ -Pi), y = qi • P\/k • P\, and | — <p is the 
azimuthal angle of the pion w.r.t. the reaction plane. They are related to the ones previously used 
by q 2 « -\Q 2 and £ « x B /(2 - x B ). 

A calculation of the matrix element squared results into 

— ,^ A9 , , = ^4 ^ b(1 ~^ { 8(1 " x B )H ud *H ud - -^-x 2 B S ud *S ud - Ax 2 B Re (n ud *S ud 
dQ 2 dx B d\A 2 \d<p Q 8 N 2 (2 - x B ) 2 \ v ; 2M 2 B B V 

-4x B VT^f^J 1 - sin <p Im j, (15) 
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where we have chosen the transverse polarization for the target proton S = (0, cos<P, sin$, 0) along 
y-axis, <P = tt/2, so that it does not possess the longitudinal component in the laboratory frame 
with z-axis along the lepton beam. The azimuthal angle <p is between the proton spin and the 
projection of the pion momentum on the (x, |/)-plane. For the minimal momentum transfer we 
use the approximation A^ in « —M 2 x\j{\ — x B ), which is valid for x-qM 2 / Q 2 <^ 1 — x-q. We also 
dropped kinematical x|M 2 /Q 2 corrections in the cross section (f[5l). The functions Ti, and £ are 
introduced as follows 



(£, A 2 ) = J o du jf_ i dx 4>MT ud (u, x, I ~ d \ [x, C, A 2 ) . (16) 




The conversion from leptoproduction to photoproduction, da£ , with longitudinally polarized 
photons e L is done by multiplication of the result (|T5| ) by kinematical factors, namely, 



+ + ( \e L -j\ 2 \ ( \L-j\ 2 d*k' \ + 1 ii x-q Q 2 

^ =d ° \-^K)\k.P x 2E^Y) = d ° —^—ydx^dQ 2 - (17) 

Here, to get the last equality we used \e L ■ j\ 2 = 4Q 2 , \L ■ j\ 2 = 16(1 — y)/y 2 , d?k'/(2E') = 
Txy/(2x-o) dx-sdQ 2 and the known definition of y resulting from the ratio of the flux factors. 



3 Next-to-leading order corrections. 

Let us turn to the evaluation of the NLO corrections to the hard-scattering amplitude. As we noted 
above in section the hard-scattering amplitude T ud can be deduced from that one computed in 
the perturbative approach to the pion form factor^ by the replacement v — > and restoration 
of the imaginary part according to the causal Feynman prescription. Thus, we decompose the 
amplitude as 



C F a s (n R )^-\Q u T ( i0, 1 u 



2£ I V 2 £ 

where 



Q, /J F , Vr) - QdT f - iO, u 



T (v, u | Q, ii F , fi R ) = i- 1 1 + ^M T (D (l7j u | Gj ^) + o ^ | . 



(19) 



Note that due to this iO-prescription, T ud satisfy an unsubtracted dispersion relation that is con- 
sistent with that one for the 7 L p — > 7i + n amplitude [|]]. Analogous to the pion form factor we 

2 One has to use those results where the symmetry properties of the distribution amplitudes, the hard-scattering 
subprocess is convoluted with, was not used. 
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extracted the LO coefficient function from the NLO amplitude T^ 1 '. In our consequent presenta- 
tion we will drop for simplicity the dimensionfull arguments, i.e. Q, the factorization scale /ip and 
the renormalization scale /xr. 

The one-loop correction for the pion form factor has been calculated by different authors in 
dimensional regularization, however, different renormalization and factorization prescriptions were 
applied || ^, [K| [II], [H| [R|. The occurring differences were clarified in Refs. JIT , [I~2|j . Indeed, 
if one takes into account errors, which are pointed out in Ref. [fl2 |, the results given in Refs. 
H [12] can be transformed to those in Refs. || [10], [11], They are given in the MS scheme and 



the subtractions of infinities are done in a way that respects the universality of the distribution 
amplitude and the running coupling in this scheme. Thus, we have 



C F T F + (3 TP+ (C F -^)T FA 



(20) 



with 



T F = [3 + \n(vu)] In f % ) + - ln 2 (t;w) + 3 ln(im) 



In w 



2 



1 ^ f \ 5 
2 lnW "6' 



2(1 



FA 



+ 



+ 



Li 2 (l - f ) - Li 2 (t>) + ln(l - v) In 
1 



u 



l-u 

{v + u — 2vu) ln(l — v) + 2tm ln(t>) 



(f — w) 5 
(1 - w)v 2 + (1 - u)u 2 



Inu 



2(1 



u 
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v — u 



[ln(l - u) ln(u) - Li 2 (l - v) + Li 2 (v)] ^ + «-> u}, (21) 



where Li 2 (w) = 
defined by Ca ~- 
lT F N f - fN c . 



/ U — ln(l — x) is the Euler dilogarithm. The colour factors are conventionally 



N c , C F 



N 2 - 
2N, 



- and the first coefficient of the QCD beta function reads /3q 



4 Leading vs. next-to-leading order observables. 

To provide estimates for the cross section and the transverse single spin asymmetry as well as 
to discuss the scale setting ambiguities we use the following models for the GPDs. For H ud we 
assume a factorizable ansatz of (x,£) and A 2 dependence valid at small A 2 . In terms of the double 
distribution F(y,z), we have 

H ud (x, & A 2 ) = G(A 2 ) ^ dy dz5(y + £z - x)F ud (y, z). (22) 

J-i J-i+\v\ 
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Here the form factor has a dipole parametrization G(A 2 ) = (1 — A 2 /m\) 2 with m 2 A = 0.9 GeV 2 
and unit normalization. The double distribution is modeled according to Ref. |l4j 



F(y, z) = {Aq(y)9(y) - Aq(-y)9(-y)} n(\y\,z), n(y, z) = 7 l „ yj ■ (23) 



3(l-y) 2 _-z 2 

yf 



Assuming an SU(2) symmetric sea, i.e. Aq ud {y) = Au(y) — Ad(y) = 0, we have Aq ud (y) = 
Au va \(y) — Ad va x(y). For the evaluations given below we use GSA forward parton densities at 
Q 2 = 4 GeV 2 ||16|| . Note that the simple ansatz (|22"D for H ud is chosen in such a way that 
constraints arising from the reduction to the forward limit and the sum rule for the lowest moment 
are satisfied. 

For E ud we adopt the pion pole-dominated ansatz 

E ud (x,£,A 2 ) = ^(A 2 )^^0.(^) , (24) 

with ^(A 2 ) taken in our estimates in the form given in Eq. (39) of Ref. In the vicinity of the 
pion pole this form factor reads -^(A 2 — > m 2 ) = 4gAM 2 / (m 2 — A 2 ). For the pion distribution 
amplitude (f) n we will take for simplicity its asymptotic formj] 

<f> v (u) = <p asy (u) = 6m(1 - u), < u < I. (25) 

In the following we will discard the evolution effects of the GPDs for simplicity since they are 
small for the models and initial scale we chose. Especially, they are negligible for the asymptotic 
pion distribution amplitude since they arise solely from two-loop (and higher) effects. Same applies 
to E since it is proportional to the asymptotic distribution amplitude, too, i.e. £E ud (x,£) oc 



4> n ((x + £)/(2£)). The neglected effects are of order 1% or so | 17fl . For H the evolution is more 
prominent since it shows up already at LO. However, it is still small for the change of scale from 
4 GeV 2 to 10 GeV 2 and result in a 5 — 8% change in the shape [[L7j . 



Let us now discuss the scale setting issues in the NLO coefficient function. A truncation of 
the perturbative series at a given order of coupling, here at NLO, causes a residual dependence on 
the factorization and renormalization scales. Obviously, this fact implies scale setting ambiguities, 
which result into uncertainties for the theoretical predictions. There are different suggestions to 
solve the scale setting problem with the aim to minimize the unknown higher order corrections. 
Unfortunately, they can provide quite different theoretical predictions for observables. 

We start with the discussion of the hf dependence, which is intimately related to the evolution 

of the distribution amplitudes/GPDs. Note that this dependence completely disappears in S, if we 

3 It results into a good agreement of the theoretical predictions for the photon-to-pion transition form factor with 
experimental data. However, one should be aware that other quite different distribution amplitudes are consistent 
with the data (see for instance fill). 



7 



assume (as we do) that 0^ coincides with the asymptotic distribution amplitude and one neglects 
its NLO evolution. The coefficient of \n(Q 2 /n F ) in Eq. (|2l| ) is merely the LO evolution kernel 



(convoluted with the LO coefficient function) and the asymptotic distribution amplitude is its 
eigenfunction with zero eigenvalue, e.g. 



dx E ud (x, (3 + 2 In = 0. 



In 7i the term proportional to In (Q 2 /fi F ) survives now in the convolution with H. However, 
as we pointed out above the evolution at LO is weak and the NLO corrections are small Wh - 
This implies that the /z^-dependence in the 7i amplitude is very weak. Therefore, we simply set 
fip = Q in what follows. A detailed discussion of other plausible settings for \i F runs beyond the 
scope of the present paper. 

Next we turn to the fiR scale. The variation of the amplitude with this scale appears at LO 
from the change of a s (fin). At NLO this scale ambiguity is expected to be smaller since this 
change is partially canceled by the variation of the term proportional to (3q in Eq. (|21|). For 
reasons explained below, we will concentrate here only on two possibilities: 

• The renormalization scale is set equal Q, fin = Q. 



The Brodsky-Lepage-Mackenzie (BLM) scale setting prescription |fl9 



For the naive setting we observe in both 7i, separately for its real and imaginary part, and 
in £ function a cancelation of two large contributions: a positive one proportional to {—(3q) and 
a negative one proportional to Cp. The term proportional to {Cp — Ca/2) in Eq. (|^) is not 
numerically important since it is relatively suppressed by 1/N 2 . Taken together they result for 
Nf = 3 into a large positive effect. 

The naive scale settings may not be quite appropriate by different reasons. The renormalization 
and factorization scales reflect quite different types of perturbative corrections Q The first one 
arises from the renormalization of ultraviolet divergences and generates the running of a s , while the 
second one comes from the factorization of collinear singularities and provides the evolution of the 
distribution amplitudes and GPDs. The former originates from the geometric series with expansion 
terms having definite sign. On the other hand, the latter is of the evolution and remnant Sudakov 
type effects ||20|| . The exponentiated Sudakov corrections have a sign alternating expansion. Thus, 
while both partially cancel at NLO, it is expected that they will amplify at NNLO. Furthermore, 
since the hard external scattering scale is partitioned among a number of parton participants in 
the hard-scattering, their resulting virtuality is much softer than the external hard scale. Thus, 



4 In the consequent discussion we follow analogous considerations for the pion form factor in Ref. 



s 



the coupling in the quark-gluon emission vertex has to be taken at a mean virtuality of e.g. hard 
gluon. A generalization of this idea results into a sensible prescription to absorb all effects coming 
from the running of the coupling, i.e. terms proportional to the /3-function, into the scale of the 
coupling. This is the BLM procedure |I9] . Since the trace anomaly of the energy momentum 



tensor is proportional to the /9-function, the resulting series in a s , if using this procedure, formally 
coincides with the perturbation theory in the conformally invariant (massless) QCD. 

Since in general we have to deal with two different amplitudes defined in Eq. ([if]), which contain 
a real and imaginary part, the consequent application of the BLM setting procedure requires the 
introduction of different scales. We separately determine the scales for the real and imaginary 
parts of the functions TL and 8 from the conditions 

. — / dx^{u Q u — P 4\ ~ Q d — cT ~„ M> A =0. (26) 

q u J—i 1 £ - x - i0 £ + x -iO 11 E ud J V 1 

Note that the scales will depend on the shape of the distribution amplitude and the GPDs. 

For 8 there is no imaginary part since E, as defined in Eq. (0), is concentrated solely in the 
exclusive domain. For its real part due to x — > —x symmetry both the Q u and Qd contributions 
in Eq. flTSP are the same and one immediately finds that the BLM scale 

/4 = Q 2 e~ u/3 (27) 



coincides with the one of the pion form factor [21, 20] for the asymptotic distribution amplitude. 



To discuss the scale setting in 7i, let us first consider the properties of H. From Eqs. 
with Aq = it is obvious that the function vanishes for x < — £, H ud (x < — £, £) = 0. From this 
we conclude that the imaginary part of Qd contribution vanishes. As we mentioned above, for H 
the scale will be different for the imaginary and real parts and will depend on the skewedness 

/4=<2 2 e->«>. (28) 

For the imaginary part one gets 

M q » ^-/■fr l-^'.fl/^tt.fl (29 ) 
6 2£ Ji Z-x 

The scale for the real part of the 7i- contribution is governed by the function 

_ 19 J-i dx {Q u ^- x In g - Q d ^- X In Jgg + g (QJ($ - x) - Q d 5({j + x))} H ud (x, Q 

(30) 

where the symbol PV stands here for the Cauchy principal value prescription. 



For the models we are using in all of the three cases e - -^ is below 0.1 in a wide interval of 
Bjorken variable 0.1 < Xb < 0.5. This means that the argument of the coupling is driven into 
the infrared region. There are experimental indications that the coupling is frozen at such a low 



scale ||22|| . Indeed, in exclusive processes, such as nucleon form factors, fixed angle proton-proton 
elastic scattering etc., the data at large scales are consistent with the predictions of dimensional 
counting rules. On the other hand, the perturbative leading twist QCD analyses coincides with 
the dimensional counting rules, however, is proportional to powers of a s , e.g. a 2 /Q A and a^/t 8 , 
respectively. Since higher order analyses favour a low scale in the coupling, we may conclude that 



the latter is a slowly varying function in this domain [32]. For our purposes we choose a frozen 
coupling with a s ([i R )/ir = 0.1 for fj, R < 1 GeV 2 ET], EDI instead of using three different soft scales 
in the analytical coupling, e.g. with gluon mass fi R — > fi 2 R + Am 2 . 

Now we are in a position to present numerical estimates of the observables at LO and NLO. 
Due to the asymptotic form of the distribution amplitude, which we use, the integration w.r.t. the 
momentum fraction u can be done analytically: 

T ud (x } = C du ^(u)T ud (u, x, (31) 



C F a s ^ R )^ \q u T (L-± - iOj - Q d T - iO 



where 



T(v) = C du(T y (u)T(v,u) (32) 
Jo 

1 Jo , a s(VR) ( n rrF ( ..\ , orr0f..\ , ( r* C A\ n - F A / 



with 



, 3 + ^ (C F r F {v) + (3 T%v) + [C F -^f) T F \v) )+0(a 



T F = - [3 + 2 In v] In ( ^] + - Inv [3 + \nv) 3 77 



2 L " ' ~\t4j ' 2~"~ L ~ ' ~"~ J 2(1 -v) 4 

T FA = -1 - 6[\nv + 2 ln(l - v)] + 12(1 - 3v) (hi 2 (v) - Li 2 (l - v) + In 



+6(1 - Qv + Qv 2 ){3U 3 (v) + 3Li 3 (l - v) - In — ^— U 2 (v) + ln 2 (l - v) \nv 

1 — v 



TV 2 



— [ln^ + 21n(l- V )]j>, (33) 

and Li 3 (t>) = Jq <fuLi 2 («)/«. 

The LO predictions were done for the running coupling with = Q, Aq° d = 220 MeV and 
Nf = 3. At NLO we use, as discussed above, two scale setting procedures: the naive \ir = Q, 
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and BLM one, with the running a s (Q 2 ) in the first and a fixed a s /ir = 0.1 below 1 GeV 2 in the 
second case, respectively. Our LO predictions given in Fig. [| (a, b) are in agreement with recent 
analyses in Refs. || [| [7|. However, they are plagued by large uncertainties from the higher order 
corrections. It is interesting to note that taking forward parton distribution as a model for the 
GPD, one gets almost identical results for the cross section. 
The transverse single spin asymmetry defined by 

ft*, daf r** daf \ { r 2 * , daf \~ l , OA . 
= U d\A2\d~tp ~ I dfA^jrf^ j U 'djA^\dip ) (34) 

is studied numerically in Fig. |] (c, d) for A 2 = —(0.1, 0.3) GeV 2 . Our leading predictions shows 
the same large asymmetry as has originally been observed in Ref . || . As our studies demonstrate 
this observable is less sensitive to higher order effects. 



5 Discussion and conclusions. 

We have discussed in the present Letter the cross section and transverse proton single spin asym- 
metry in the hard exclusive production of pions. For the longitudinal polarization of the virtual 
photon the amplitude is dominated by twist-two contributions, i.e. lowest Fock components in 
the pion and p — > n transition amplitudes. We evaluated the NLO corrections to the hard par- 
ton subprocess basing on the available result for the pion form factor and estimated their size in 
physical observables. We found that the ambiguity in the renormalization scale setting (together 
with the freezing of the coupling in the infrared region) results into large uncer tainties in the 
theoretical predictions, (T^lo — Tlo)/Tlo, of the amplitude. After application of the BLM scale 
setting prescription, we observed a large reduction in the magnitude of the cross section related to 
the effects of Sudakov double logarithms. A deeper insight into the structure of these corrections 
is highly desirable. 

The theoretical uncertainty in the factorization procedure on the amplitude level is translated 
into large variations of the physical cross section. However, we found that the single spin asymme- 
try, given by the ratio of the Fourier coefficients of the cross section, is a 'good' observable since 
the ambiguities due to the truncation of the perturbative series approximately cancel. Thus, the 
perturbative predictions for this quantity are rather stable. The NLO effects result into cor- 
rections to the LO prediction for 0.1 < xb < 0.5. For the assumed models of the non-perturbative 
functions the asymmetry is big, being as large as 60%, and is sensitive to the pion pole dominated 
GPD E. Thus, in view of its advantages being rather insensitive to the higher order corrections, it 
turns out to be an appropriate quantity for experimental studies at Jefferson Lab and HERMES. 
The large NLO corrections may indicate that the application of the perturbative approach to the 
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Figure 4: The leading twist predictions for the unpolarized photoproduction cross section 
d&l + /d\ A 2 1 at Q = 10 GeV 2 for the GPD models specified in Eqs. (H,H are shown for A 2 = A^ in 
and A 2 = —0.3 GeV 2 in the panels (a) and (b), respectively. In (c) and (d) we display the trans- 
verse proton single spin asymmetry A± for A 2 = —0.1 GeV 2 and A 2 = —0.3 GeV 2 , respectively. 
The solid, dashed and dash-dotted curves represent the LO and NLO with the naive and BLM 
scale setting, respectively. 

meson production cross section is legitimate at a rather large momentum transfer. However, for 
the asymmetry, which depends at leading twist only logarithmically on Q 2 , one expects its earlier 
validity due to an essential cancelation of higher order perturbative and, hopefully, also power 
suppressed corrections, due to their intrinsic interrelation in field theoretical treatment, recall 
renormalons in this respect. 

Let us stress that the longitudinal proton single spin asymmetry, measured at HERMES ||, 
arises at the twist-three level and requires a separate analysis. The result we have presented here 
are given for the transversely polarized proton target. The NLO analysis addressed in our study 
can be extended to a large set of meson leptoproduction amplitudes with the quark dominated 
short-distance subprocess, e.g. £p — > £'n°p, £p — > £'p + n, etc. 

We would like to thank J.C. Collins, P. Kroll, J. Smith, M. Stratmann for useful conversations 
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J. Smith for the hospitality extended to him at the C.N. Yang Institute for Theoretical Physics 
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